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We observe that, within the effective generating function formalism for the implementation of 
canonical transformations within wave mechanics, non-trivial canonical transformations which leave 
invariant the form of the Hamilton function of the classical analogue of a quantum system manifest 
themselves in an integral equation for its stationary state eigenfunctions. We restrict ourselves 
to that subclass of these dynamical symmetries for which the corresponding effective generating 
functions are necessarily free of quantum corrections. We demonstrate that infinite families of such 
OS ' transformations exist for a variety of familiar conservative systems of one degree of freedom. We show 

OS ' how the geometry of the canonical transformations and the symmetry of the effective generating 

function can be exploited to pin down the precise form of the integral equations for stationary state 
eigenfunctions. We recover several integral equations found in the literature on standard special 
functions of mathematical physics. We end with a brief discussion (relevant to string theory) of the 
generalization to scalar field theories in 1+1 dimensions. 
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> ! I. INTRODUCTION 

<N 
O 

Given the formal similarities between quantum mechanics and the Hamiltonian formulation of classical mechanics, 
l/^ I it is not surprizing that there have been several attempts to define within quantum mechanics transformations anal- 
ogous to the canonical transformations underpinning the powerful Hamilton- Jacobi method. The founding fathers of 
On ■ quantum mechanics were content to identify the quantum analogues of canonical transformations as unitary transfor- 
mations of the position and momentum operators which preserve the canonical commutation relations fjLES. However, 
to enhance the scope of applications, various operator-based extensions of this notion of a quantum canonical trans- 
Qh formation have been proposed ||,|J . The path- integral formulation of quantum mechanics (and quantum field theory) 
' with its c-number representation of co-ordinates and momenta has also prompted more ambitious constructions [gj , 
but these are not without their complications Q (because of the ambiguous interrelationships between the momenta 
Pi and co-ordinates qi introduced in the discretization of the path integrals Jf/j). Finally, there are the realizations of 
canonical transformations within quantum theory as integral transforms |^,^||T 

It has been observed j^] that, at least for time- independent problems, there is a simple implementation of canonical 
transformations within wave mechanics via an integral transform of the form 



Mq) = n a j eW F ^$ a (Q)dQ. (1.1) 

Here, <j> a {q) and Q a (Q) denote the wavefunctions of the same state \a) for two choices of generalized co-ordinate q 
and Q, n a is a state-dependent relative normalization (independent of q and Q) and, to leading order in h, the state- 
independent quantum (or effective) generating function F(q, Q) coincides with the generating function T(q, Q) of the 
canonical transformation q,p — ► Q, P within classical mechanics. As we demonstrate below (in Sec. |0|), this relation 
between classical and quantum generating functions respects the algebra of composition of canonical transformations 
and, under certain broad conditions (spelt out in Sec. |l|), the quantum corrections to T{q, Q) vanish. (Unless stated 
otherwise, we assume throughout this wo rk that we are dealing with systems of one degree of freedom but use of the 
quantum canonical transform in Eq. (1.1) is not limited to such systems |iT^ . |lT| , p^| , ^3| .) 



Quantum realizations of canonical transformations provide elegant solutions to eigenvalue problems in quantum 
mechanics |[3 14 ,lj|. They also offer a platform for the exploration of a novel kind of issue, namely the impact at 



the quantum level of symmetries which hold in the full phase space at the classical level. A topical example would 
be the investigation of duality in field theories beyond the classical level |ll],[l6|]. In this paper, we consider the 
implications of the existence of non-trivial canonical transformations q,p — ► Q, P in the full phase space such that 
the transformed Hamiltonian function (or Kamiltonian) JC(Q,P) = TL(q(Q , P) , p(Q , P)) is the same as the original 
Hamiltonian function H., i.e. IC(Q,P) — H.(Q.P). To this end, we adopt the effective generating function formalism 
of§. 
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Canonical transformations which preserve the form of a Hamilton function TL have a dynamical significance within 
the classical mechanics of systems of one degree of freedom: they map trajectories of such systems (level curves of 
H.) onto themselves and, thus, amount to evolutions in time. Within wave mechanics, we have for a form-preserving 
transformation q,p — > Q, P the equality <f> a (Q) = <j> a {Q) (because the Hamiltonian operators for q and Q are identical 



in form) and the integral transform in Eq. (1.1) reduces to an integral equation for stationary state eigenfunctions. 

The point of departure for the present study is the conjecture that the class of (classical) generating functions 
T{q, Q) of form-preserving canonical transformations may include members for which quantum corrections vanish. 
We restrict ourselves to Hamiltonians describing a particle of mass m in the potential V(q). We demonstrate (in 
Sec. |n] below) that, in addition to the free theory (V = 0), there are six distinct choices of potential V(q) for which 
infinite families {F^q^Q)} of non-trivial form-preserving and correction- free generating functions exist (jj, labels the 
members of these families). These potentials include some of the most ubiquitous (notably the linear and quadratic 
potentials) and the integral equations implied 

= JV„(/*) I e (W " (9 ' Q Va(Q)dQ (1.2) 



apply to several standard special functions of mathematical physics: in the notation of |17[| , the parabolic cylinder 
functions D n (integer order), the Airy function Ai, the modified Bessel functions K a (imaginary order), the Mathicu 

functions ce r and se r and the modified Mathieu functions Mc?' and Msi 1 ^. 

The major preoccupation of this paper is to show that not only can the kernel e^ i ' n ' :F ' t ^ q '^ in these integral equa- 
tions be obtained via consideration of canonical transformations, but so can the reciprocals j\f a (ll) of its eigenvalues 
(modulo, in some cases, a phase). Two distinct and complementary methods apply. For those potentials for which the 
canonical transformations form an abelian group (the linear and quadratic potentials), the composition of canonical 
transformations implies a functional relation for j\f a (p) which determines it up to a phase. For the other potentials, we 
can take advantage of a remarkable symmetry in the dependence of the corresponding generating functions J-p,(q, Q) 
on fx, q and Q (described in Sec. |v|). 



We close in Sec. VI with some remarks on lines of investigation suggested by this work. 



II. THE QUANTUM CANONICAL TRANSFORM 

Let {4> a (q)} denote the complete set of stationary wavefunctions of a quantum system with Hamiltonian operator 
h (q, and let H ^Q, j-§q^ be the realization of the Hamiltonian of the system for another choice of generalized 

co-ordinate Q and {$ Q (Q)} the corresponding complete set of stationary wavefunctions. The quantum cano nica l 



transform introduced in || relates the Q (g)'s to the <& Q ((3)'s via an integral relationship of the form in Eq. ([Lip . 
The quantum generating function F(q, Q) is fixed by the requirement that 4> a {q) and <fr a (Q) are eigenfunctions of the 
same complete set of commuting observables with exactly the same set of quantum numbers a. The state-dependent 
relative normalization n a is chosen so that the normalizations of 4> a (q) and <fr a (Q) are compatible. 

For a system with a non-degenerate energy spectrum (like those considered in Sees. IV and |v| below), the restriction 



on F(q, Q) reduces to the condition that cf> a (q) and $ a (Q) are eigenfunctions of h and H, respectively, with the same 



energy E a . Substituting for cf> a in h<j) a — E a <p a using Eq. (1.1) and then replacing the product E a $? a by H$? a , wc find 



that, after the requisite number of integration by parts and appealing to the completeness of the ^q's, this condition 
implies F(q, Q) should satisfy 

provided the endpoint terms generated in the integration by parts vanish. These terms take the form of the bilinear 
combination 

p {i/h)F{q,Q) _ (J( ^_ p (i/n)F{q,Q) ( ~ o) 

dQ dQ 1 ' 

for the Hamiltonian 

n 9 \ ^ d 2 



H { Q ^)=-^w +v{Q) (2 - 3) 

considered below. For bound states, the vanishing of the wavefunction and its first derivative at infinity guarantees 



that the bilinear concomitant in Eq. (2.2) is zero 
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The relation of the quantum generating function to a generating function f c \(q, Q) of a canonical transformation 
within the Hamiltonian formulation of classical mechanics can be brought out by adopting for F(q, Q) an expansion 
in powers of ih: 



F(q,Q) = £V n fa,Q)(ift)' 



n=0 



If we take H (q, j§q^ to be given by Eq. ( 2^ ) and h (q, ||^J to be given by 



2m dq 2 



then substitution of this expansion into Eq. (2.1) yields for !Fq 

( m 



i 

2m V 



v(q) 



1 

2m 



and for the other J^n's (n > 0) 



E( dT k dT n - k _ dT k dT n - k 
\ dq dq dQ dQ 



dfo\ 
dQ) 



V(Q) 



dQ 2 



(2.4) 



(2.5) 



Equation ( |2.4| ) is automatically satisfied if we identify To(q,Q) as the classical generating function of a canonical 
transformation (q,p) — > (Q,P) for which the original Hamiltonian function TC(q,p) = p 2 /(2m) + v(q) and the 
transformed Hamiltonian function (or Kamiltonian) /C(Q,P) = P 2 /(2m) + V(Q). If d 2 J r /dq 2 = d 2 J r /dQ 2 , i.e. T 
is of the form 



fo(q,Q) = J r +(q+)+J r -(q-), 



(2.6) 



where q± = (q ± Q)/2 and J-+ and T- are arbitrary functions, then Eq. (2.5) implies that the quantum corrections 
•7"n ( n ^ 1) can be taken to be zero. 

The parallel between quantum and classical generating functions also extends to the composition of transformations. 
Let Fi(q, q^) and -F^fe, Q) denote the quantum generating functions for the canonical transformations (q,p) — ► (q%,Pi) 
and (qi,Pi) — ► (Q,P), respectively. The exact relation among these generating functions and the quantum generating 
function F c (q,Q) for the composition (q,p) — ► {Q,P) reads 



n (1) n (2) 



^l/h)[F 1 (q,q i )+F 2 (q u Q)] 



(c)M/h)F c (q,Q) 



(2.7) 



To identify the relation between the classical contribution f c to F c and the classical con tributions to F\ and F% 
(assumed to be /i and fz, respectively), we can evaluate the integration over qi in Eq. ( |2.7| ) in the stationary phase 
approximation. Retaining only the terms most singular in H 7 we obtain 



n (,c) e (i/h)f (q,Q) 



2nh 



!„(l)„(2) p (i/?i)/sO?,Q) 



(2.8) 



where 



with qi chosen so that 



fs(q, Q) = fi{q, qi) + fz{qi, Q) 



7r (h(q,qi)+f2(q i ,Q)) 
oqi 



= 



(2.9) 



(2.10) 



Qi=qi 



and we have made the generically valid assumption that 



d 2 h , d 2 h 
dq 2 dq 2 
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is non-zero. (For simplicity, we have also assumed that there is only one stationary point.) Interpretation of Eq. (2.S 



is complicated by the unknown, but, in general, singular d epe ndence of the relative normalizations ria \ nffl and tv* 
on h. Nevertheless, the essential singularities in h in Eq. ( |2.8| ) can only match if the dependence on q and Q cancels 
— i.e. f c and f s differ at most by a constant. Exactly the same relation between f c and f s is implied by the canonical 
formalism of classical mechanics. We have that 



dq dq V % dqi J i _ dq dq 
dQ dQ \ d qi dqj, _ dQ dQ 

1 1i=<li 



= +P, 



= -P, 



i.e. the partial derivatives of f s coincide with those of f c . 

Further parallels with generating functions in the canonical formalism of classical mechanics are discussed in Gq,13l 



III. INVARIANT HAMILTONIAN FUNCTIONS 



Quantum generating functions F(q, Q) which give rise to integral equations can be obtained by determining classical 
generating functions Fo{q, Q) of the form in Eq. ( |2.6| ) which induce canonical transformations q,p — > Q, P such that 
the Kamiltonian function fC(Q,P) = H(q(Q, P),p(Q, P)) is the same as the original Hamiltonian function Tt : i.e. 
K.{Q,P) = H(Q,P). (Below, we drop the subscript 0, denoting T by T.) 

We consider a canonically conjugate pair of variables q and p for which the Hamiltonian function 



Inspection of the relations for the momenta in terms of a classical generating function J-(q, Q) 



(3.1) 



shows that at least for the case of the free theory (V = 0) there are non-trivial canonical transformations under which 
the form of the Hamiltonian is unchanged: for the generating functions .F f * c = J--(qJ) and F f ~ cc — {q+), where 
J- + and JF_ are arbitrary, the transformed generalized momentum P = +p and P = —p. respectively, from which the 
form invariance of the free Hamiltonian Hf ree = p 2 /(2m) is obvious. 
More generally, substituting for the momenta in 

v 2 P 2 
^+V(q) = — + V(Q) 
2m Am 



using Eq. (|3.l| ), we deduce that the following relation must hold between V and T±: 

I 



2m 



T'{q+)T'_{q-) - V(q + - q_) - V(q + + q_). 



(3.2) 



To proceed, we assume that T'_ and V are analytic and we expand both sides of Eq. ( |3.2| ) in powers of g_ to obtain 
(* = 0,1,2,...) 



-^^;(x)^i fc+1) (o) = \(-\f - i]vW(x), 

lm 



(3.3) 



where we have set q+ = x. If we ignore the possibility that V'(x) = (in which case we recover the results given 
above for the free theory), then Eq. ( |3.3| ) for k = 1 implies that \x = !F_ (0) ^ and 

4m. 



-V'(x). 



(3.4) 



Substituting Eq. (3.4) into Eq. (3.3). it reduces to the simultaneous requirements that odd derivatives of V are given 
by 
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(3.5) 



and that T— is e ven (so that the derivatives T_ (0) = for n odd). 
From Eq. (3.4), we immediately have that 



4m 
p 



where we have dropped an irrelevant constant of integration. The implications of Eq. (| 
whether or not the third derivative V^ 3 -* vanishes. 



3) for V and T- depend on 



If V^(x) ^ 0, then Eq. (|J) implies that p = (0) / ' T {T) (0) ^ 0, V^{x) = pV^(x) and g {2k) (0) = p k ^p 
for k > 1. Thus, since /? may be of either sign, the potential can either be the combination of hyperbolic functions 
( P = +/3 2 >0) 

V+ (x) = A cosh (3x + B sinh fix 
or the combination of sinusoidal functions (p = — 1 < 0) 

V_ (x) — A cos (3x + B sin /3a;, 
where /?, A and £? are arbitrary constants. The corresponding forms of T- are (up to arbitrary additive constants) 

(x) — cosh [3x 

P 



and 



TI(x) 



if 



cos (3x, 



respectively. If we assume that V^(x) = 0, then Eq. ( |3.5| ) implies that, with the exception of T_ (0), all the even 

derivatives ^ r i 2n ' ) (0) = 0. Discarding arbitrary additive constants, the potential is of the quadratic form V(x) — 
Ax 2 + Bx, where, as above, A and B are arbitrary constants, and is th e qu adratic T-(x) = p,x 2 /2. 

Not only are there several classes of potential V(q) compatible with Eq. ( |3.2| ), but also, for each class of potentials, 
there is an infinite family of canonical transformations distinguished by different values of the parameter p, which is 
not fixed by the considerations above. By translation of the origin or the inversion q — > —q or translation followed by 
inversion, all members of the classes of non-trivial potentials identified above can be reduced to one of those in Table |. 
We also include the families of generating functions of canonical transformations which leave the corresponding 
Hamiltonians unchanged. In what follows, we take A > 0, a choice which embraces the physically more interesting 
scenarios. 

Below [in Sec. @], we shall have cause to invoke the limit p — > oo. From Eqs. ( |3.l| ) and ( |3.4[ ), the difference in 
momenta 



P-P=— V(q + ) 



(3.6) 



The difference in the coordinates Q — q (= — 2g_) in terms of p + P can be obtained by inversion of the relation 



P + p = TL{q-) 



(3.7) 



implied by Eq. (3.1). Together, Eqs. (3.6) and (3.7) imply that, for our choices of J--, both P — p = and Q — q are 



of orde r l/ p for large p. [In the case of the sinusoidal potential, use of the principal value of arcsin in the inversion 
of Eq. (3.7) is understood.] Hence, in the limit p — > oo, all the canonical transformations of interest reduce to the 
identity transformation. 
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IV. INTEGRAL EQUATIONS FOR THE QUADRATIC AND LINEAR POTENTIALS 



It is natural to ask whether any of the families of canonical transformations we have identified are groups? At 
the level of generating functions, this question translates into whether when two generating functions and 
of a particular family in Table | are combined according to Eq. (2.£), the resulting generating function is also 

a member of that family (to within an additive constant). This is not automatic since thes e fa milies of generating 
functions have been constructed by restricting their members to be of the special form in Eq. (2.6) and this constraint 
cannot, in general, be respected by the prescription in Eq. (2.E). However, it turns out that, because of the simple 
geometrical operations in phase space effected by members of the families of canonical transformations corresponding 
to the quadratic and linear poten tial s, these two families do form groups. We are able to exploit the abelian character 
of these groups to convert Eq. (2J5) into a functional relationship for the reciprocals Af a (n) of eigenvalues and so 
deduce the integral equations for the quadratic and linear potentials. 



A. The quadratic potential 

The canonical transformations which leave the Hamiltonian function associated with this potential unchanged are 
linear transformations of the phase plane and so must be as either a rotation or an area-preserving shear or an area- 
preserving squeeze [jl9). The observation that, for the rescaled canonical variables q s = (mA) 1 / 4 ^ and p s = p/ (mA) 1 / 4 , 
the Hamiltonian function 



suggests that it should be possible to decompose these transformations in terms of a rotation as 



Q 
P 



_ ( (mA)- 1 / 4 

(mA) 1 / 4 



cos o sm t 
— sin 9 cos ( 



(mA) 1 / 4 


° ^ 


1 





(mA)" 1 / 4 J 


P 



where the angle of rotation 9 is dependent on the choice [i. In fact, this does prove to be the case with the free 
parameter fi uniquely related to the angle of rotation 9 by 



a = —2V\m cot -. 
p 2 

The totality of these linear transformations thus constitutes a faithful matrix representation of the rotation group 
SO(2). 

It is convenient to work with generating functions parametrized by 9 instead of namely 



F(q, Q\9) = ^mcu [2 esc 9qQ - cot 9{q 2 + Q 2 )] , 



(4.1) 



where lo = \J A/m. We write the corresponding integral equation for eigenfunctions {tp n (q)} (n a non-negative integer) 
of the harmonic oscillator hamiltonian 



2m dq 2 



2 2 

—muoq 
2 



as 



(4.2) 



Since h\ lo commutes with the parity operator P and has a non-degenerate spectrum, its eigenfunctions are auto- 
matically either even or odd. Although the generating function F(q,Q\9) has been constructed only with a view to 
ensuring that the integral on the righthand-side of Eq. (^^) is an eigenfunction of /iho with the same eigenenergy as 
ijj n (q), the dependence of F{q, Q\9) on q and Q guarantees that this integral has also the same parity as tp n (q). The 
generating function pertinent to the sinusoidal potential shares this property. 

Use of 9 in Eq. (4.2) simplifies the treatment of the composition of transformations. We can immediately identify 
the quantum generating function F c (q,Q) for the composition of two transformations with generating functions 



G 



F\ = F(q,qi\0i) and F 2 — F(qi,Q\02), respectively, as F c (q,Q) = F(q,Q\9 1 +#2)- Accordingly, in this context, 
Eq. reads 

iV n (0i)iV n (0 2 ) />/»>F(«.«l*>+^«.0lfe)]d© = 7V„(0i + ^ 2 ) e (W(9,Q|e 1 +e 2 ) ) 
which, on evaluation of the gaussian integral over qi, becomes 



N n (0 1+ 2 ) = XI^J^^NMN^). (4.3) 



The functional relation in Eq. (4.3) has the solution 



where the coefficient c n is arbitrary. 

The coefficient c n can be fixed by appealing to the fact that, under the inversion q,p — ► —q, —p [corresponding to 
the choice of = it in F(q, Q\0)], the eigenfunctions ip n (q) transform in a well-defined manner: ip n {~q) = (— l) n i/) n {Q)- 
From Eqs. (Q and Q» 

lim Ar n (0)e (i/a)F(9 ' Q|e) = ie^fa + Q), 

which, on substitution in Eq. ( |4.2|) , implies ip n (q) = ie Cri7V ^p n (—q). Consistency with the parity properties of the 
ip n (q)'s is achieved by taking c„ = — (n + l/2)i. 

In its final form, our integral equation for the eigenfunctions ipn(l) of the harmonic oscillator reads 



To make the connection with known results, we note that we may read off from Eq. (|4.5| ) that the expansion coefficients 
of jV ((9)e( l /' i - )F ( 9 ' ( 3l e ) in the orthonormal basis {ip n (Q)} are e m6 ip n (q) and, hence, construct the identity 

^/^cscSgQ-cot^+Q 2 )] = [^1^0 y e +i(n+l/2)0 ^ ( ^ (Q) _ 

V mwi 

n 

Invoking the relation of the tp n (q)'s to the parabolic cylinder function D n (x) fl7| l, namely 

^n(?) = -^=D n (y / 2muj/hq), 



we recover one of the addition theorems for the -D ra 's given in chapter II of p0|| . (To obtain this addition theorem in 
the precise form given in p0|, we must set e ie — itan0, yj mui / hq — e 47r / 4 A and yj muj / TiQ = e l37r / 4 /i.) 



B. The linear potential 

The relevant canonical transformations involve an area-preserving shear of the phase plane coupled with a shift of 
the origin: 



P 



1 -2v 
1 



1m\v 



I 



where we have introduced the parameter v = 2/ 1 \x which makes the algebra of these transformations more transparent. 
(The factor of 2 in the definition of v is a matter of convenience.) Because of the shift in origin, it is perhaps not 
obvious that this class of geometrical operations should form a group, but, in fact, the composition 



T T — T 



(4.6) 
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so that the transformations T v are a representation of an abelian affine group. 

Let F v (q, Q) denote the generating function of these canonical transformations when v and not \x is adopted as the 
free parameter, i.e. 

Fu(q,Q) =-m\v{q + Q) + ±-(q-Q)\ 



and let i/je{q) (E real) denote the eigenfunction of 

^linear ; . 

Zm oq 



n 2 d 2 , 



with eigenenergy E. Identifying F l5 F 2 and F c in Eq. fl2.7| ) with F Ul , F V2 and F V1+V2 , respectively, we find that the 
reciprocals Ne(v) of the eigenvalues in the integral equation involving the eigenfunction iPe{i) must satisfy 

N E {vi + v 2 ) = V^iJ^^e-^^^^+^NEMNEi^). (4.7) 



Equation (4.7) determines Ne(v) to be of the form 

N E {V) = r—— L 

V VKTllV 

The E dependence of the coefficient ce can be pinned down by using the relation of eigenfunctions of non-zero energy 
ipE(q) to the zero energy eigenfunction ipQ (q) : 

Mq)=VeMq~E/X), (4.8) 

where t\e is a constant of modulus unity which we assume below is absorbed into the definition of iIje(q) with an 
appropriate choice of the phase. 

Using Eq. (4.8) to substitute for iPe(Q) and ipE(q) in the integral equation 



^ E (q) = N E {v) J e^ F ^ Q ^j E (Q)dQ, (4.9) 
we find, after the change of variable Q — > Q' = Q — E / A, 

Mq ~ E/X) = N E {v)e- amEv ' h [ e^ F ^- E / X 'Q^MQ')dQ' = ^e^/^^ - E/X), 



where to obtain the last equality we have invoked Eq. (4.9) again. The choice ce = 2mE is indicated. 

An independent check of these results is given by working with the momentum space equivalent of Eq. (4.9): 

$b(p) = N E {v) J K(p, P)$ E (P)dP, (4.10) 

where the kernel 

K(p, P)= ( e -to/n e m)M*Q) e WQ/nMtoL = ^m) e -^'^ n ^ +p ^ 8(P -p- 2mXv) 
J 2irn 

and the Fourier transform 

MP) = -^= / e- iw/S fe(?)d9 = e-* E P/^Mp). 

The delta function in K (p, P) enfo rces the relation between the momenta in the canonical transformation generated 
by F v (q, Q) and reduces Eq. ( 4.10 ) to the algebraic equation 

Mp)e- tp3/(6mXh) = N E (v)V4^e~^ 2mE »~ m2x2 » 3 /VM P + 2mXv)e~^ +2mXv ?l^ mXh \ 
since ip (p) = Ce +lpa /( 6mXh \ where C is a normalization constant, we recover the results above for Ne(v)- 
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The zero energy eigenfunction ^0(3) [and hence all the other eigenfunctions iPe(q)] is related to the Airy function 
Ai(x) @: 

where 7 = (2mA/?! 2 ) 1 / 3 . [For the sake of definiteness, the norm aliz ation is fixed so that ('0-eIV'e) = — E').] Thus, 
in terms of the Airy function Ai(x), our integral equation Eq. (4.9) amounts to the relation 

Ai(z) = -jL=e- is "I 12 / e ^(*+x)/2+(*-x)V(4s)] Ai(x)dX; 
where we have set x — 75, X — 7Q and s = h^ 2 v. 



V. INTEGRAL EQUATIONS FOR THE OTHER POTENTIALS 



In Sec. Ill we found that, for potentials V(x) for which V^ 3 ^ 7^ 0, non-trivial form-preserving correction-free gener- 
ating functions J-^(q, Q) only exist if V"(x) — pV(x) (p a constant). A related implication is that the dependence of 
these generating functions on q and Q must be such that 



d 2 Tn r d 2 T, 



dq 2 



dQ 2 



(5.1) 



The dependence on p is such that 



which, setting p = p{z) = p$e^P z , becomes 



dp 



pT, 



(5.2) 



The similarity of Eqs. (5.1) and (5.2) suggests that it should be possible to treat z in T^u^q, Q) in formally the same 
way as the generalized co-ordinates q and Q. In fact, we find that, with appropriate choices of po (which are listed in 
Table ||, 



•?>(*) Q) = •?>(</) (2, Q) = -7>(Q) z ) 

confirming that the roles of z and q (or 2 and Q) may be interch ang ed. 

Used in conjunction with the integral equation Eq. (1.2), Eq. (5.3) implies that 

M a (p(z)) 

Thus, the reciprocals Af a (p) of eigenvalues are given to within a constant C' a by 



(5.3) 



= / e (VS)^ ( *)(9.Q)^ a (Q)dQ = / e «/»)^(.)(*.«)^ Q (Q)dQ = -^M 



(/*(«))' 



M a {p) 



where 2 = = hi(p,/ p.0) 1 / The constant C Q can be fixed by the requirement that 

lim M a {p)e {l/h):F ^ Q) = <% - Q), 



(5.4) 



(5.5) 



reflecting the fact that, in the limit p — > 00, we recover the identity transformation from T^q^Q) (cf. the end of 
Sec. 

To establish the limit on the lefthand-side of Eq. (5.5), we consider the integral 

3 (</ft}*>(9.«J)$(Q)dQ, (5.6) 

where $((3) is a suitable test function, and apply the method of stationary phase to extract the leading contribution 
to the integral when p^5> 1. 

By way of illustration, we now discuss the cases of the exponential and sinusoidal potentials in more detail. 
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A. The exponential potential 



The hamiltonian 



Exponential- 2 m flg 2 + 2 



has eigenf unctions of energy E k = h 2 k 2 /(2m) given up to a normalization constant by 

ipk(q) cx K l(k/a) (y^e a y[ha 2 ]^j , (5.7) 



where K a denotes a modified Bessel function of order a (section 9.6 in JL7|). Substitution of Eq. (5.7) into Eq. (5.4) 
yields for the reciprocals of eigenvalues 

A4(/i) = C k /K Kk/a) (^/[Atha 2 }) , (5.8) 

where the constant C k has still to be determined by consideration of the fi — ► oo limit. 

Asymptotic analysis for //> 1 via the method of stationary phase implies that, to leading order, the integral 



J V A 1 



Thus, Eq. (|5.5|) is satisfied provided 

Hm J^e-r'^KV) = " (5.9) 

which, on the substitution of Eq. (5.S) and use of the leading term in the asymptotic expansion of K a (x) for x 3> 1, 
reduces to the requirement that C k — a/2. 

In terms of the variables y = V m\ae aq / (ha 2 ), Y = v m\ae a ® / (ha 2 ), p — ik/a and w = /i / (Ahia 2 ) , our integral 
equation reads 

oo 

2K p (w)K p (y) = J e -^ w +^y/ Y + Y /y^K p (Y)^, 
o 

which coincides formally with Eq. 6.653.2 in pl| (after the change of integration variable Y — ► x = wy/Y). 



B. The sinusoidal potential 

We confine our attention to the denumerable set of eigenfunctions ip s (q) (s — 0, ±1, ±2, . . .) of 

7 h 2 d 2 A 

/^sinusoidal = ~ + T 1 ? cos ^ a Q 

2m oq z 4or 

which are related to the Mathieu functions ce r and se r : 

!C s ce s (aq,S) s = 0,1,2,... 
(5.10) 
C s se\ s \(aq, 6) s = -1, -2, . . . , 

where the dimensionless rescaled strength of the potential S = mA/(4?i 2 a 4 ) and C s denotes a normalization constant. 
The choice of the range of integration in the integral equation for the ?/> s 's 

A(q) = N s ([i) [ e^o>Q^ s (Q)dQ 
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is dictated by the consideration that the bilinear concomitant in Eq. (2.2) vanishes. This can be achieved by exploiting 
the periodicity of the generating function T^(q, Q) (period 2n/a in q or Q) and the eigenfunctions ip s (q) (period ir/a 
for s even and 2tt / 'a for s odd). Thus, we take the range of integration to be over one period of the generating function 
from to +27r/a. (Use of these non-symmetric limits facilitates the comparison with the integral equations tabulated 
in ch. 20 of 



Substituting for the wavefunction in Eq. (5.4) using Eq. ( 5.1C ) and invoking the proportionality of ce r (— ix, rj) and 
se r (—ix,rj) to the modified Mathieu functions Mcv • (x,rj) and M Sr v)i respectively, we find that 



M s (1) (\n[^/V4m\},6 



(5.11) 



where C s is independent of p, and m] 1 -* (x, ?y) denotes the modified Mathieu function Meg 1 ' (x, rj) for s > and the 
modified Mathieu function Ms^) (x,tj) for s < 0. 

In the asymptotic analysis of the integral in Eq. ( |5.6| ) for large /x, we encounter in the present case two points of 
stationary phase: Q = Qi = q + 0(p~ 2 ) and Q = Q2 = q — a q -Kja + 0(/i~ 2 ), where o q denotes the sign of q — w/a. 
To leading order, the asymptotic expansion reads 



2-K/a 



(i/h)^(g,Q) 



$(g)dg 



I8irh 
A 6 



(5.12) 



To proceed, it is necessary to recognise that the choice of the appropriate space of test functions $(q) depends on the 
properties of the eigenfunctions {tjj s (q)} under consideration. Accordingly, $((7) is drawn from either of two spaces 
S p (p = 0, 1): a space So of periodic functions of period Tr/a appropriate to the eigenfunctions ip s for even s, and: a 
space S\ of periodic functions of period 2n/a appropriate to the eigenfunctions ip s for odd s. Furthermore, paralleling 
the property that ipsil ~ a q^lo) = {~ ^ sW) 1 we must require th at te st functions drawn from S\ are such that 
$(<7 — <j q iT/a) = — $(<?). Thus, for test functions drawn from S p , Eq. ( 5.12 ) becomes 



27r/a 



$(Q)dQ ~ 2i~P* 



cos[p/{4ha 2 ) - (p + 1/2)tt/2]$(£ ? ) 



(5.13) 



so that, for a suitable choice of N s (fi), lim p >oc N s {p)e^ l / h ^ :F ^ q ^ can have the sifting property expected. 

To leading order, the asymptotic expansion of Afj (mj/i/V 4mA] , 5) in the limit of large fj, (> 0) is [from the real 
part of Eq. (20.9.1) in 0] 



m\], 5 



I8ha 2 



lS [ M /(4aa 2 )-(p+l/2)7r/2], 



(5.14) 



where p = (1) for s even (odd). [Despite appearances, the righthand-side of Eq. (5.14) is real-valued consistent with 
the reality of Afj 1 ^ (x, rf) for real-valued arguments x and 77.] Combining Eqs. (5.11), (5.13) and ( |5.14 ) in Eq. fl5.5| ), 
we conclude that C s = i' s ' a/ (27r) . 

Introducing £ = ln[///v4mA] and the variables u = q/a and U = Q/a, the integral equation for the ip s 's reads 

271 



i> s (u/a) 



i\s\ 



1 



2^M S (1) (C,5) 



? — 2iV^(cosh £ cos u cos t/ +sinh £ sin u sin J7) 



^ s (C//a)df/, 



which for £ real is tantamount to the complex conjugate of Eqs. (20.7.34) and (20.7.35) in pj 
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VI. DISCUSSION 



In this paper, we have explored the implications for quantum mechanical systems of the existence of non-trivial 
canonical transformations which leave the form of the Hamilton function of the corresponding classical system invari- 
ant. Such dynamic al s ymmetries of a classical system manifest themselves in a linear homogeneous integral equation 
of the kind in Eq. ( |l.2| ) for the stationary state eigenfunctions of the quantum mechanical system. We have seen that 
the quantum canonical transform introduced in |9[| is ideally suited to the purpose of tying down the features of these 
integral equations. 

We have not exhausted the full range of integral equations which can be constructed by invoking the quantum 
canonical transform. For the families of canonical transformations which we have considered which do not constitute 
Abelian groups, there is a (possibly infinite) sequence of integral equations corresponding to the repeated composition 
of these transformations. Consideration of the composition of a transformation to angle-action variables with its 
inverse could give rise to still more integral equations. However, the most interesting line of further investigation in 
our opinion is the extension of the work in this paper to quantum field theory. 

There are some almost immediate parallels of our results for theories of a scalar field </?(ct, t) in 1+1 dimensions (a 
denotes the spatial dimension and t the time in natural units such that Ti = 1 = c). This is in part a consequence 
of the fact that the (first quantized) Hamiltonian functionals H[ip,n] for these theories are not too dissimilar in form 
from the Hamilton functions Ti. we have considered: 



H[<p,n] = \l \tt 2 + {dip/ d<jf \ dff + y V(ip)da 



where tt is the field momentum conjugate to <p and the "potential density" V describes the self-coupling of ip. If we 
consider canonical transformations ip, tt — > <I>, n induced by generating functionals of the form [tp± = (ip ± $)/2] 

F[<p,*] = J ^da + J [F + (<p+)+F-(<p-)]dcr, 

then we find that a sufficient condition for the transformed Hamilton functional to be of the same form as the original 
Hamilton functional is that 

^^ = V( V+ - t p-)-V( kV+ + „.) (6.1) 
provided F± either vanish or are periodic at the endpoints of the integration over a. Apart from an inessential factor 



of m, Eq. (p.l| ) is formally identical to the condition established in Sec. HI for the invariance of Hamilton functions 
under canonical transformations [Eq. ( |3.2| )] with V(ip) replacing the potential V(q) and p±{(J,t) the combinations 
q± = (q ± Q)/2. In the Schrodinger representation |22| of the corresponding (second) quantized field theories, we 
have, putting aside the issue of renormalization, a class of integral equations for the wavefunctionals of the form 

where, in terms of the generating functions listed in Table |, 

f 5$ f 
F M b,*] = / ^da + J ^,$)da 

and M a is a functional of fi which is now a function of a and r. 



[1] 
[2] 
[3] 
[4] 

[5] 



M. Born, W. Heisenberg and P. Jordan, Z. Phys. 35, 557 (1926). 
P. A. M. Dirac, Proc. R. Soc. London Ser. A 110, 561 (192 6). 
A. Anderson, Ann. Phys. 232, 292 (1994), |hep-th/93050 54 



H. W. Lee and W. S. l'Yi, Phys. Rev. A 51, 982 (1994), [hep-th/9406106 



A. K. Kapoor, Phys. Rev. D 29, 2339 (1984); ibid. 30, 1750 (1984). 



12 



[6] 
[7] 
[8 
[9 
[10 

[11 
[12 
[13 

[11 
[15 
[16 

[ir 

[18 

[19 
[20 
[21 

[22 



M. Blau, E. Keski-Vakkuri and A. Niemi, Phys. Lett. B 246, 92 (1990). 
C. Garrod, Rev. Mod. Phys. 38, 483 (1966). 
J. Hietarinta, Phys. Rev. D 25, 2103 (1982). 
G. I. Ghandour, Phys. Rev. D 35, 1289 (1987). 

T. L. Curtright and C. Zachos, Phys. Rev D 49, 5408 (19 94), |hep-th/940l"006 . 
Y. Lozano, Phys. Lett. B 364, 19 (1995), |hep-th/9503045j 
A. J. Bordner, J. Math. Phys. 38, 3427 (1997). 

K. F. Elshamandi and G. I. Ghandour, The R adial Generators for Central Potentials, submitted for p ublication. 

A. Anderson, J. Math. Phys. 35, 6018 (19 94), |hep-th/9310168j ; ibid. 37, 1349 (1996), |hep-th/950709^ . 
W. S. 1'Yi, Phys. Rev. A 53, 1251 (1996), ^iep-th/9508'027t 

K. Sfetsos, Phys. Lett. B 432, 365 (1998), hep-th/98030T9| . 

Handbook of Mathematical Functions, eds. M. Abramowitz and I. A. Stegun (NBS, Washington, 1964). 

S. M. Mady, On the Effective Generating Functions for Quantum Canonical Transformation, M.Sc. thesis (Kuwait University, 

1988), unpublished. 

I. Percival and D. Richards, Introduction to Dynamics (Cambridge University Press, Cambridge, 1982), ch. 2. 
P. M. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw-Hill, New York, 1953). 

I. S. Gradshteyn and I. M. Rhyzik, Tables of Integrals, Series and Products, 5th ed. (Academic Press, San Diego, 1994). 

B. Hatfield, Quantum Field Theory of Point Particles and Strings (Addison- Wesley, Redwood City, 1992), ch. 10. 



TABLE I. 


Standard forms of potentials V and the corresponding 


(classical) generating functions Tp, 




V(g) 




Quadratic 




-t (q + Q)2 + -s {q - Q)2 


Sinusoidal 


— cos 2aq 

4a 2 


mX , u , „, 
— cos a(q + Q ) — cos a(q — Q ) 


Even Hyperbolic 


— ^- cosh 2aq 
4a 2 


— -^4 cosn a(q + Q) + cosh a(q — Q) 
[ia A Aa J 


Linear 


Xq 


(q + Q) + s(q-Q) 

fl o 


Exponential 


— e 2aq 
2a 


2mX ea{q+Q)+ u, 
fia 4a 2 


Odd Hyperbolic 


— sinh 2aq 
2a H 


— ^H!^. sinh a(q + Q) + cosn a (<7 — Q) 
fj,a 4ar 



TABLE II. Parameters of the transformation z(u.) = p ln(/x//xo) 



V(g) Mo 0> 

Sinusoidal 2V mA ia 

Even hyperbolic 2 V mXi a 

Exponential 4-V mXai a 

Odd hyperbolic 2\/2mXai a 
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